\(\int \frac {1}{x^2 (a+b x^4+c x^8)} \, dx\) [325]

Optimal result

Integrand size = 18, antiderivative size = 363 \[ \int \frac {1}{x^2 \left (a+b x^4+c x^8\right )} \, dx=-\frac {1}{a x}-\frac {\sqrt [4]{c} \left (1-\frac {b}{\sqrt {b^2-4 a c}}\right ) \arctan \left (\frac {\sqrt [4]{2} \sqrt [4]{c} x}{\sqrt [4]{-b-\sqrt {b^2-4 a c}}}\right )}{2\ 2^{3/4} a \sqrt [4]{-b-\sqrt {b^2-4 a c}}}-\frac {\sqrt [4]{c} \left (1+\frac {b}{\sqrt {b^2-4 a c}}\right ) \arctan \left (\frac {\sqrt [4]{2} \sqrt [4]{c} x}{\sqrt [4]{-b+\sqrt {b^2-4 a c}}}\right )}{2\ 2^{3/4} a \sqrt [4]{-b+\sqrt {b^2-4 a c}}}+\frac {\sqrt [4]{c} \left (1-\frac {b}{\sqrt {b^2-4 a c}}\right ) \text {arctanh}\left (\frac {\sqrt [4]{2} \sqrt [4]{c} x}{\sqrt [4]{-b-\sqrt {b^2-4 a c}}}\right )}{2\ 2^{3/4} a \sqrt [4]{-b-\sqrt {b^2-4 a c}}}+\frac {\sqrt [4]{c} \left (1+\frac {b}{\sqrt {b^2-4 a c}}\right ) \text {arctanh}\left (\frac {\sqrt [4]{2} \sqrt [4]{c} x}{\sqrt [4]{-b+\sqrt {b^2-4 a c}}}\right )}{2\ 2^{3/4} a \sqrt [4]{-b+\sqrt {b^2-4 a c}}} \]

[Out]

-1/a/x-1/4*c^(1/4)*arctan(2^(1/4)*c^(1/4)*x/(-b-(-4*a*c+b^2)^(1/2))^(1/4))*(1-b/(-4*a*c+b^2)^(1/2))*2^(1/4)/a/
(-b-(-4*a*c+b^2)^(1/2))^(1/4)+1/4*c^(1/4)*arctanh(2^(1/4)*c^(1/4)*x/(-b-(-4*a*c+b^2)^(1/2))^(1/4))*(1-b/(-4*a*
c+b^2)^(1/2))*2^(1/4)/a/(-b-(-4*a*c+b^2)^(1/2))^(1/4)-1/4*c^(1/4)*arctan(2^(1/4)*c^(1/4)*x/(-b+(-4*a*c+b^2)^(1
/2))^(1/4))*(1+b/(-4*a*c+b^2)^(1/2))*2^(1/4)/a/(-b+(-4*a*c+b^2)^(1/2))^(1/4)+1/4*c^(1/4)*arctanh(2^(1/4)*c^(1/
4)*x/(-b+(-4*a*c+b^2)^(1/2))^(1/4))*(1+b/(-4*a*c+b^2)^(1/2))*2^(1/4)/a/(-b+(-4*a*c+b^2)^(1/2))^(1/4)

Rubi [A] (verified)

Time = 0.25 (sec) , antiderivative size = 363, normalized size of antiderivative = 1.00, number of steps used = 8, number of rules used = 5, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.278, Rules used = {1382, 1524, 304, 211, 214} \[ \int \frac {1}{x^2 \left (a+b x^4+c x^8\right )} \, dx=-\frac {\sqrt [4]{c} \left (1-\frac {b}{\sqrt {b^2-4 a c}}\right ) \arctan \left (\frac {\sqrt [4]{2} \sqrt [4]{c} x}{\sqrt [4]{-\sqrt {b^2-4 a c}-b}}\right )}{2\ 2^{3/4} a \sqrt [4]{-\sqrt {b^2-4 a c}-b}}-\frac {\sqrt [4]{c} \left (\frac {b}{\sqrt {b^2-4 a c}}+1\right ) \arctan \left (\frac {\sqrt [4]{2} \sqrt [4]{c} x}{\sqrt [4]{\sqrt {b^2-4 a c}-b}}\right )}{2\ 2^{3/4} a \sqrt [4]{\sqrt {b^2-4 a c}-b}}+\frac {\sqrt [4]{c} \left (1-\frac {b}{\sqrt {b^2-4 a c}}\right ) \text {arctanh}\left (\frac {\sqrt [4]{2} \sqrt [4]{c} x}{\sqrt [4]{-\sqrt {b^2-4 a c}-b}}\right )}{2\ 2^{3/4} a \sqrt [4]{-\sqrt {b^2-4 a c}-b}}+\frac {\sqrt [4]{c} \left (\frac {b}{\sqrt {b^2-4 a c}}+1\right ) \text {arctanh}\left (\frac {\sqrt [4]{2} \sqrt [4]{c} x}{\sqrt [4]{\sqrt {b^2-4 a c}-b}}\right )}{2\ 2^{3/4} a \sqrt [4]{\sqrt {b^2-4 a c}-b}}-\frac {1}{a x} \]

[In]

Int[1/(x^2*(a + b*x^4 + c*x^8)),x]

[Out]

-(1/(a*x)) - (c^(1/4)*(1 - b/Sqrt[b^2 - 4*a*c])*ArcTan[(2^(1/4)*c^(1/4)*x)/(-b - Sqrt[b^2 - 4*a*c])^(1/4)])/(2
*2^(3/4)*a*(-b - Sqrt[b^2 - 4*a*c])^(1/4)) - (c^(1/4)*(1 + b/Sqrt[b^2 - 4*a*c])*ArcTan[(2^(1/4)*c^(1/4)*x)/(-b
+ Sqrt[b^2 - 4*a*c])^(1/4)])/(2*2^(3/4)*a*(-b + Sqrt[b^2 - 4*a*c])^(1/4)) + (c^(1/4)*(1 - b/Sqrt[b^2 - 4*a*c]
)*ArcTanh[(2^(1/4)*c^(1/4)*x)/(-b - Sqrt[b^2 - 4*a*c])^(1/4)])/(2*2^(3/4)*a*(-b - Sqrt[b^2 - 4*a*c])^(1/4)) +
(c^(1/4)*(1 + b/Sqrt[b^2 - 4*a*c])*ArcTanh[(2^(1/4)*c^(1/4)*x)/(-b + Sqrt[b^2 - 4*a*c])^(1/4)])/(2*2^(3/4)*a*(
-b + Sqrt[b^2 - 4*a*c])^(1/4))

Rule 211

Int[((a_) + (b_.)*(x_)^2)^(-1), x_Symbol] :> Simp[(Rt[a/b, 2]/a)*ArcTan[x/Rt[a/b, 2]], x] /; FreeQ[{a, b}, x]
&& PosQ[a/b]

Rule 214

Int[((a_) + (b_.)*(x_)^2)^(-1), x_Symbol] :> Simp[(Rt[-a/b, 2]/a)*ArcTanh[x/Rt[-a/b, 2]], x] /; FreeQ[{a, b},
x] && NegQ[a/b]

Rule 304

Int[(x_)^2/((a_) + (b_.)*(x_)^4), x_Symbol] :> With[{r = Numerator[Rt[-a/b, 2]], s = Denominator[Rt[-a/b, 2]]}
, Dist[s/(2*b), Int[1/(r + s*x^2), x], x] - Dist[s/(2*b), Int[1/(r - s*x^2), x], x]] /; FreeQ[{a, b}, x] &&  !
GtQ[a/b, 0]

Rule 1382

Int[((d_.)*(x_))^(m_)*((a_) + (c_.)*(x_)^(n2_.) + (b_.)*(x_)^(n_))^(p_), x_Symbol] :> Simp[(d*x)^(m + 1)*((a +
b*x^n + c*x^(2*n))^(p + 1)/(a*d*(m + 1))), x] - Dist[1/(a*d^n*(m + 1)), Int[(d*x)^(m + n)*(b*(m + n*(p + 1) +
1) + c*(m + 2*n*(p + 1) + 1)*x^n)*(a + b*x^n + c*x^(2*n))^p, x], x] /; FreeQ[{a, b, c, d, p}, x] && EqQ[n2, 2
*n] && NeQ[b^2 - 4*a*c, 0] && IGtQ[n, 0] && LtQ[m, -1] && IntegerQ[p]

Rule 1524

Int[(((f_.)*(x_))^(m_.)*((d_) + (e_.)*(x_)^(n_)))/((a_) + (b_.)*(x_)^(n_) + (c_.)*(x_)^(n2_)), x_Symbol] :> Wi
th[{q = Rt[b^2 - 4*a*c, 2]}, Dist[e/2 + (2*c*d - b*e)/(2*q), Int[(f*x)^m/(b/2 - q/2 + c*x^n), x], x] + Dist[e/
2 - (2*c*d - b*e)/(2*q), Int[(f*x)^m/(b/2 + q/2 + c*x^n), x], x]] /; FreeQ[{a, b, c, d, e, f, m}, x] && EqQ[n2
, 2*n] && NeQ[b^2 - 4*a*c, 0] && IGtQ[n, 0]

Rubi steps \begin{align*} \text {integral}& = -\frac {1}{a x}+\frac {\int \frac {x^2 \left (-b-c x^4\right )}{a+b x^4+c x^8} \, dx}{a} \\ & = -\frac {1}{a x}-\frac {\left (c \left (1-\frac {b}{\sqrt {b^2-4 a c}}\right )\right ) \int \frac {x^2}{\frac {b}{2}+\frac {1}{2} \sqrt {b^2-4 a c}+c x^4} \, dx}{2 a}-\frac {\left (c \left (1+\frac {b}{\sqrt {b^2-4 a c}}\right )\right ) \int \frac {x^2}{\frac {b}{2}-\frac {1}{2} \sqrt {b^2-4 a c}+c x^4} \, dx}{2 a} \\ & = -\frac {1}{a x}+\frac {\left (\sqrt {c} \left (1-\frac {b}{\sqrt {b^2-4 a c}}\right )\right ) \int \frac {1}{\sqrt {-b-\sqrt {b^2-4 a c}}-\sqrt {2} \sqrt {c} x^2} \, dx}{2 \sqrt {2} a}-\frac {\left (\sqrt {c} \left (1-\frac {b}{\sqrt {b^2-4 a c}}\right )\right ) \int \frac {1}{\sqrt {-b-\sqrt {b^2-4 a c}}+\sqrt {2} \sqrt {c} x^2} \, dx}{2 \sqrt {2} a}+\frac {\left (\sqrt {c} \left (1+\frac {b}{\sqrt {b^2-4 a c}}\right )\right ) \int \frac {1}{\sqrt {-b+\sqrt {b^2-4 a c}}-\sqrt {2} \sqrt {c} x^2} \, dx}{2 \sqrt {2} a}-\frac {\left (\sqrt {c} \left (1+\frac {b}{\sqrt {b^2-4 a c}}\right )\right ) \int \frac {1}{\sqrt {-b+\sqrt {b^2-4 a c}}+\sqrt {2} \sqrt {c} x^2} \, dx}{2 \sqrt {2} a} \\ & = -\frac {1}{a x}-\frac {\sqrt [4]{c} \left (1-\frac {b}{\sqrt {b^2-4 a c}}\right ) \tan ^{-1}\left (\frac {\sqrt [4]{2} \sqrt [4]{c} x}{\sqrt [4]{-b-\sqrt {b^2-4 a c}}}\right )}{2\ 2^{3/4} a \sqrt [4]{-b-\sqrt {b^2-4 a c}}}-\frac {\sqrt [4]{c} \left (1+\frac {b}{\sqrt {b^2-4 a c}}\right ) \tan ^{-1}\left (\frac {\sqrt [4]{2} \sqrt [4]{c} x}{\sqrt [4]{-b+\sqrt {b^2-4 a c}}}\right )}{2\ 2^{3/4} a \sqrt [4]{-b+\sqrt {b^2-4 a c}}}+\frac {\sqrt [4]{c} \left (1-\frac {b}{\sqrt {b^2-4 a c}}\right ) \tanh ^{-1}\left (\frac {\sqrt [4]{2} \sqrt [4]{c} x}{\sqrt [4]{-b-\sqrt {b^2-4 a c}}}\right )}{2\ 2^{3/4} a \sqrt [4]{-b-\sqrt {b^2-4 a c}}}+\frac {\sqrt [4]{c} \left (1+\frac {b}{\sqrt {b^2-4 a c}}\right ) \tanh ^{-1}\left (\frac {\sqrt [4]{2} \sqrt [4]{c} x}{\sqrt [4]{-b+\sqrt {b^2-4 a c}}}\right )}{2\ 2^{3/4} a \sqrt [4]{-b+\sqrt {b^2-4 a c}}} \\ \end{align*}

Mathematica [C] (verified)

Result contains higher order function than in optimal. Order 9 vs. order 3 in optimal.

Time = 0.02 (sec) , antiderivative size = 71, normalized size of antiderivative = 0.20 \[ \int \frac {1}{x^2 \left (a+b x^4+c x^8\right )} \, dx=-\frac {1}{a x}-\frac {\text {RootSum}\left [a+b \text {$\#$1}^4+c \text {$\#$1}^8\&,\frac {b \log (x-\text {$\#$1})+c \log (x-\text {$\#$1}) \text {$\#$1}^4}{b \text {$\#$1}+2 c \text {$\#$1}^5}\&\right ]}{4 a} \]

[In]

Integrate[1/(x^2*(a + b*x^4 + c*x^8)),x]

[Out]

-(1/(a*x)) - RootSum[a + b*#1^4 + c*#1^8 & , (b*Log[x - #1] + c*Log[x - #1]*#1^4)/(b*#1 + 2*c*#1^5) & ]/(4*a)

Maple [C] (verified)

Result contains higher order function than in optimal. Order 9 vs. order 3.

Time = 0.09 (sec) , antiderivative size = 63, normalized size of antiderivative = 0.17

method	result	size
default	\(-\frac {\munderset {\textit {\_R} =\operatorname {RootOf}\left (c \,\textit {\_Z}^{8}+\textit {\_Z}^{4} b +a \right )}{\sum }\frac {\left (\textit {\_R}^{6} c +b \,\textit {\_R}^{2}\right ) \ln \left (x -\textit {\_R} \right )}{2 \textit {\_R}^{7} c +\textit {\_R}^{3} b}}{4 a}-\frac {1}{a x}\)	\(63\)
risch	\(-\frac {1}{a x}+\frac {\left (\munderset {\textit {\_R} =\operatorname {RootOf}\left (\left (256 a^{9} c^{4}-256 b^{2} c^{3} a^{8}+96 b^{4} c^{2} a^{7}-16 b^{6} c \,a^{6}+b^{8} a^{5}\right ) \textit {\_Z}^{8}+\left (80 a^{4} b \,c^{4}-120 a^{3} b^{3} c^{3}+61 a^{2} b^{5} c^{2}-13 c \,b^{7} a +b^{9}\right ) \textit {\_Z}^{4}+c^{5}\right )}{\sum }\textit {\_R} \ln \left (\left (\left (1152 a^{9} c^{4}-1184 b^{2} c^{3} a^{8}+456 b^{4} c^{2} a^{7}-78 b^{6} c \,a^{6}+5 b^{8} a^{5}\right ) \textit {\_R}^{8}+\left (332 a^{4} b \,c^{4}-487 a^{3} b^{3} c^{3}+245 a^{2} b^{5} c^{2}-52 c \,b^{7} a +4 b^{9}\right ) \textit {\_R}^{4}+4 c^{5}\right ) x +\left (64 a^{8} c^{4}-112 a^{7} b^{2} c^{3}+60 a^{6} b^{4} c^{2}-13 a^{5} b^{6} c +a^{4} b^{8}\right ) \textit {\_R}^{7}+\left (4 a^{3} b \,c^{4}-a^{2} b^{3} c^{3}\right ) \textit {\_R}^{3}\right )\right )}{4}\)	\(305\)

[In]

int(1/x^2/(c*x^8+b*x^4+a),x,method=_RETURNVERBOSE)

[Out]

-1/4/a*sum((_R^6*c+_R^2*b)/(2*_R^7*c+_R^3*b)*ln(x-_R),_R=RootOf(_Z^8*c+_Z^4*b+a))-1/a/x

Fricas [B] (verification not implemented)

Leaf count of result is larger than twice the leaf count of optimal. 5758 vs. \(2 (281) = 562\).

Time = 0.53 (sec) , antiderivative size = 5758, normalized size of antiderivative = 15.86 \[ \int \frac {1}{x^2 \left (a+b x^4+c x^8\right )} \, dx=\text {Too large to display} \]

[In]

integrate(1/x^2/(c*x^8+b*x^4+a),x, algorithm="fricas")

[Out]

Too large to include

Sympy [F(-1)]

Timed out. \[ \int \frac {1}{x^2 \left (a+b x^4+c x^8\right )} \, dx=\text {Timed out} \]

[In]

integrate(1/x**2/(c*x**8+b*x**4+a),x)

[Out]

Timed out

Maxima [F]

\[ \int \frac {1}{x^2 \left (a+b x^4+c x^8\right )} \, dx=\int { \frac {1}{{\left (c x^{8} + b x^{4} + a\right )} x^{2}} \,d x } \]

[In]

integrate(1/x^2/(c*x^8+b*x^4+a),x, algorithm="maxima")

[Out]

-integrate((c*x^6 + b*x^2)/(c*x^8 + b*x^4 + a), x)/a - 1/(a*x)

Giac [F]

\[ \int \frac {1}{x^2 \left (a+b x^4+c x^8\right )} \, dx=\int { \frac {1}{{\left (c x^{8} + b x^{4} + a\right )} x^{2}} \,d x } \]

[In]

integrate(1/x^2/(c*x^8+b*x^4+a),x, algorithm="giac")

[Out]

integrate(1/((c*x^8 + b*x^4 + a)*x^2), x)

Mupad [B] (verification not implemented)

Time = 9.01 (sec) , antiderivative size = 10509, normalized size of antiderivative = 28.95 \[ \int \frac {1}{x^2 \left (a+b x^4+c x^8\right )} \, dx=\text {Too large to display} \]

[In]

int(1/(x^2*(a + b*x^4 + c*x^8)),x)

[Out]

2*atan((((-(b^9 + b^4*(-(4*a*c - b^2)^5)^(1/2) + 80*a^4*b*c^4 + 61*a^2*b^5*c^2 - 120*a^3*b^3*c^3 + a^2*c^2*(-(
4*a*c - b^2)^5)^(1/2) - 13*a*b^7*c - 3*a*b^2*c*(-(4*a*c - b^2)^5)^(1/2))/(512*(a^5*b^8 + 256*a^9*c^4 - 16*a^6*
b^6*c + 96*a^7*b^4*c^2 - 256*a^8*b^2*c^3)))^(3/4)*(4096*a^15*c^8 - x*(-(b^9 + b^4*(-(4*a*c - b^2)^5)^(1/2) + 8
0*a^4*b*c^4 + 61*a^2*b^5*c^2 - 120*a^3*b^3*c^3 + a^2*c^2*(-(4*a*c - b^2)^5)^(1/2) - 13*a*b^7*c - 3*a*b^2*c*(-(
4*a*c - b^2)^5)^(1/2))/(512*(a^5*b^8 + 256*a^9*c^4 - 16*a^6*b^6*c + 96*a^7*b^4*c^2 - 256*a^8*b^2*c^3)))^(1/4)*
(32768*a^16*c^8 + 1024*a^12*b^8*c^4 - 12288*a^13*b^6*c^5 + 51200*a^14*b^4*c^6 - 81920*a^15*b^2*c^7)*1i + 256*a
^11*b^8*c^4 - 2816*a^12*b^6*c^5 + 10496*a^13*b^4*c^6 - 14336*a^14*b^2*c^7)*1i + 4*a^11*b*c^8*x)*(-(b^9 + b^4*(
-(4*a*c - b^2)^5)^(1/2) + 80*a^4*b*c^4 + 61*a^2*b^5*c^2 - 120*a^3*b^3*c^3 + a^2*c^2*(-(4*a*c - b^2)^5)^(1/2) -
13*a*b^7*c - 3*a*b^2*c*(-(4*a*c - b^2)^5)^(1/2))/(512*(a^5*b^8 + 256*a^9*c^4 - 16*a^6*b^6*c + 96*a^7*b^4*c^2
- 256*a^8*b^2*c^3)))^(1/4) - ((-(b^9 + b^4*(-(4*a*c - b^2)^5)^(1/2) + 80*a^4*b*c^4 + 61*a^2*b^5*c^2 - 120*a^3*
b^3*c^3 + a^2*c^2*(-(4*a*c - b^2)^5)^(1/2) - 13*a*b^7*c - 3*a*b^2*c*(-(4*a*c - b^2)^5)^(1/2))/(512*(a^5*b^8 +
256*a^9*c^4 - 16*a^6*b^6*c + 96*a^7*b^4*c^2 - 256*a^8*b^2*c^3)))^(3/4)*(4096*a^15*c^8 + x*(-(b^9 + b^4*(-(4*a*
c - b^2)^5)^(1/2) + 80*a^4*b*c^4 + 61*a^2*b^5*c^2 - 120*a^3*b^3*c^3 + a^2*c^2*(-(4*a*c - b^2)^5)^(1/2) - 13*a*
b^7*c - 3*a*b^2*c*(-(4*a*c - b^2)^5)^(1/2))/(512*(a^5*b^8 + 256*a^9*c^4 - 16*a^6*b^6*c + 96*a^7*b^4*c^2 - 256*
a^8*b^2*c^3)))^(1/4)*(32768*a^16*c^8 + 1024*a^12*b^8*c^4 - 12288*a^13*b^6*c^5 + 51200*a^14*b^4*c^6 - 81920*a^1
5*b^2*c^7)*1i + 256*a^11*b^8*c^4 - 2816*a^12*b^6*c^5 + 10496*a^13*b^4*c^6 - 14336*a^14*b^2*c^7)*1i - 4*a^11*b*
c^8*x)*(-(b^9 + b^4*(-(4*a*c - b^2)^5)^(1/2) + 80*a^4*b*c^4 + 61*a^2*b^5*c^2 - 120*a^3*b^3*c^3 + a^2*c^2*(-(4*
a*c - b^2)^5)^(1/2) - 13*a*b^7*c - 3*a*b^2*c*(-(4*a*c - b^2)^5)^(1/2))/(512*(a^5*b^8 + 256*a^9*c^4 - 16*a^6*b^
6*c + 96*a^7*b^4*c^2 - 256*a^8*b^2*c^3)))^(1/4))/(((-(b^9 + b^4*(-(4*a*c - b^2)^5)^(1/2) + 80*a^4*b*c^4 + 61*a
^2*b^5*c^2 - 120*a^3*b^3*c^3 + a^2*c^2*(-(4*a*c - b^2)^5)^(1/2) - 13*a*b^7*c - 3*a*b^2*c*(-(4*a*c - b^2)^5)^(1
/2))/(512*(a^5*b^8 + 256*a^9*c^4 - 16*a^6*b^6*c + 96*a^7*b^4*c^2 - 256*a^8*b^2*c^3)))^(3/4)*(4096*a^15*c^8 - x
*(-(b^9 + b^4*(-(4*a*c - b^2)^5)^(1/2) + 80*a^4*b*c^4 + 61*a^2*b^5*c^2 - 120*a^3*b^3*c^3 + a^2*c^2*(-(4*a*c -
b^2)^5)^(1/2) - 13*a*b^7*c - 3*a*b^2*c*(-(4*a*c - b^2)^5)^(1/2))/(512*(a^5*b^8 + 256*a^9*c^4 - 16*a^6*b^6*c +
96*a^7*b^4*c^2 - 256*a^8*b^2*c^3)))^(1/4)*(32768*a^16*c^8 + 1024*a^12*b^8*c^4 - 12288*a^13*b^6*c^5 + 51200*a^1
4*b^4*c^6 - 81920*a^15*b^2*c^7)*1i + 256*a^11*b^8*c^4 - 2816*a^12*b^6*c^5 + 10496*a^13*b^4*c^6 - 14336*a^14*b^
2*c^7)*1i + 4*a^11*b*c^8*x)*(-(b^9 + b^4*(-(4*a*c - b^2)^5)^(1/2) + 80*a^4*b*c^4 + 61*a^2*b^5*c^2 - 120*a^3*b^
3*c^3 + a^2*c^2*(-(4*a*c - b^2)^5)^(1/2) - 13*a*b^7*c - 3*a*b^2*c*(-(4*a*c - b^2)^5)^(1/2))/(512*(a^5*b^8 + 25
6*a^9*c^4 - 16*a^6*b^6*c + 96*a^7*b^4*c^2 - 256*a^8*b^2*c^3)))^(1/4)*1i + ((-(b^9 + b^4*(-(4*a*c - b^2)^5)^(1/
2) + 80*a^4*b*c^4 + 61*a^2*b^5*c^2 - 120*a^3*b^3*c^3 + a^2*c^2*(-(4*a*c - b^2)^5)^(1/2) - 13*a*b^7*c - 3*a*b^2
*c*(-(4*a*c - b^2)^5)^(1/2))/(512*(a^5*b^8 + 256*a^9*c^4 - 16*a^6*b^6*c + 96*a^7*b^4*c^2 - 256*a^8*b^2*c^3)))^
(3/4)*(4096*a^15*c^8 + x*(-(b^9 + b^4*(-(4*a*c - b^2)^5)^(1/2) + 80*a^4*b*c^4 + 61*a^2*b^5*c^2 - 120*a^3*b^3*c
^3 + a^2*c^2*(-(4*a*c - b^2)^5)^(1/2) - 13*a*b^7*c - 3*a*b^2*c*(-(4*a*c - b^2)^5)^(1/2))/(512*(a^5*b^8 + 256*a
^9*c^4 - 16*a^6*b^6*c + 96*a^7*b^4*c^2 - 256*a^8*b^2*c^3)))^(1/4)*(32768*a^16*c^8 + 1024*a^12*b^8*c^4 - 12288*
a^13*b^6*c^5 + 51200*a^14*b^4*c^6 - 81920*a^15*b^2*c^7)*1i + 256*a^11*b^8*c^4 - 2816*a^12*b^6*c^5 + 10496*a^13
*b^4*c^6 - 14336*a^14*b^2*c^7)*1i - 4*a^11*b*c^8*x)*(-(b^9 + b^4*(-(4*a*c - b^2)^5)^(1/2) + 80*a^4*b*c^4 + 61*
a^2*b^5*c^2 - 120*a^3*b^3*c^3 + a^2*c^2*(-(4*a*c - b^2)^5)^(1/2) - 13*a*b^7*c - 3*a*b^2*c*(-(4*a*c - b^2)^5)^(
1/2))/(512*(a^5*b^8 + 256*a^9*c^4 - 16*a^6*b^6*c + 96*a^7*b^4*c^2 - 256*a^8*b^2*c^3)))^(1/4)*1i))*(-(b^9 + b^4
*(-(4*a*c - b^2)^5)^(1/2) + 80*a^4*b*c^4 + 61*a^2*b^5*c^2 - 120*a^3*b^3*c^3 + a^2*c^2*(-(4*a*c - b^2)^5)^(1/2)
- 13*a*b^7*c - 3*a*b^2*c*(-(4*a*c - b^2)^5)^(1/2))/(512*(a^5*b^8 + 256*a^9*c^4 - 16*a^6*b^6*c + 96*a^7*b^4*c^
2 - 256*a^8*b^2*c^3)))^(1/4) - atan((((-(b^9 - b^4*(-(4*a*c - b^2)^5)^(1/2) + 80*a^4*b*c^4 + 61*a^2*b^5*c^2 -
120*a^3*b^3*c^3 - a^2*c^2*(-(4*a*c - b^2)^5)^(1/2) - 13*a*b^7*c + 3*a*b^2*c*(-(4*a*c - b^2)^5)^(1/2))/(512*(a^
5*b^8 + 256*a^9*c^4 - 16*a^6*b^6*c + 96*a^7*b^4*c^2 - 256*a^8*b^2*c^3)))^(3/4)*(4096*a^15*c^8 + x*(-(b^9 - b^4
*(-(4*a*c - b^2)^5)^(1/2) + 80*a^4*b*c^4 + 61*a^2*b^5*c^2 - 120*a^3*b^3*c^3 - a^2*c^2*(-(4*a*c - b^2)^5)^(1/2)
- 13*a*b^7*c + 3*a*b^2*c*(-(4*a*c - b^2)^5)^(1/2))/(512*(a^5*b^8 + 256*a^9*c^4 - 16*a^6*b^6*c + 96*a^7*b^4*c^
2 - 256*a^8*b^2*c^3)))^(1/4)*(32768*a^16*c^8 + 1024*a^12*b^8*c^4 - 12288*a^13*b^6*c^5 + 51200*a^14*b^4*c^6 - 8
1920*a^15*b^2*c^7) + 256*a^11*b^8*c^4 - 2816*a^12*b^6*c^5 + 10496*a^13*b^4*c^6 - 14336*a^14*b^2*c^7) + 4*a^11*
b*c^8*x)*(-(b^9 - b^4*(-(4*a*c - b^2)^5)^(1/2) + 80*a^4*b*c^4 + 61*a^2*b^5*c^2 - 120*a^3*b^3*c^3 - a^2*c^2*(-(
4*a*c - b^2)^5)^(1/2) - 13*a*b^7*c + 3*a*b^2*c*(-(4*a*c - b^2)^5)^(1/2))/(512*(a^5*b^8 + 256*a^9*c^4 - 16*a^6*
b^6*c + 96*a^7*b^4*c^2 - 256*a^8*b^2*c^3)))^(1/4)*1i - ((-(b^9 - b^4*(-(4*a*c - b^2)^5)^(1/2) + 80*a^4*b*c^4 +
61*a^2*b^5*c^2 - 120*a^3*b^3*c^3 - a^2*c^2*(-(4*a*c - b^2)^5)^(1/2) - 13*a*b^7*c + 3*a*b^2*c*(-(4*a*c - b^2)^
5)^(1/2))/(512*(a^5*b^8 + 256*a^9*c^4 - 16*a^6*b^6*c + 96*a^7*b^4*c^2 - 256*a^8*b^2*c^3)))^(3/4)*(4096*a^15*c^
8 - x*(-(b^9 - b^4*(-(4*a*c - b^2)^5)^(1/2) + 80*a^4*b*c^4 + 61*a^2*b^5*c^2 - 120*a^3*b^3*c^3 - a^2*c^2*(-(4*a
*c - b^2)^5)^(1/2) - 13*a*b^7*c + 3*a*b^2*c*(-(4*a*c - b^2)^5)^(1/2))/(512*(a^5*b^8 + 256*a^9*c^4 - 16*a^6*b^6
*c + 96*a^7*b^4*c^2 - 256*a^8*b^2*c^3)))^(1/4)*(32768*a^16*c^8 + 1024*a^12*b^8*c^4 - 12288*a^13*b^6*c^5 + 5120
0*a^14*b^4*c^6 - 81920*a^15*b^2*c^7) + 256*a^11*b^8*c^4 - 2816*a^12*b^6*c^5 + 10496*a^13*b^4*c^6 - 14336*a^14*
b^2*c^7) - 4*a^11*b*c^8*x)*(-(b^9 - b^4*(-(4*a*c - b^2)^5)^(1/2) + 80*a^4*b*c^4 + 61*a^2*b^5*c^2 - 120*a^3*b^3
*c^3 - a^2*c^2*(-(4*a*c - b^2)^5)^(1/2) - 13*a*b^7*c + 3*a*b^2*c*(-(4*a*c - b^2)^5)^(1/2))/(512*(a^5*b^8 + 256
*a^9*c^4 - 16*a^6*b^6*c + 96*a^7*b^4*c^2 - 256*a^8*b^2*c^3)))^(1/4)*1i)/(((-(b^9 - b^4*(-(4*a*c - b^2)^5)^(1/2
) + 80*a^4*b*c^4 + 61*a^2*b^5*c^2 - 120*a^3*b^3*c^3 - a^2*c^2*(-(4*a*c - b^2)^5)^(1/2) - 13*a*b^7*c + 3*a*b^2*
c*(-(4*a*c - b^2)^5)^(1/2))/(512*(a^5*b^8 + 256*a^9*c^4 - 16*a^6*b^6*c + 96*a^7*b^4*c^2 - 256*a^8*b^2*c^3)))^(
3/4)*(4096*a^15*c^8 + x*(-(b^9 - b^4*(-(4*a*c - b^2)^5)^(1/2) + 80*a^4*b*c^4 + 61*a^2*b^5*c^2 - 120*a^3*b^3*c^
3 - a^2*c^2*(-(4*a*c - b^2)^5)^(1/2) - 13*a*b^7*c + 3*a*b^2*c*(-(4*a*c - b^2)^5)^(1/2))/(512*(a^5*b^8 + 256*a^
9*c^4 - 16*a^6*b^6*c + 96*a^7*b^4*c^2 - 256*a^8*b^2*c^3)))^(1/4)*(32768*a^16*c^8 + 1024*a^12*b^8*c^4 - 12288*a
^13*b^6*c^5 + 51200*a^14*b^4*c^6 - 81920*a^15*b^2*c^7) + 256*a^11*b^8*c^4 - 2816*a^12*b^6*c^5 + 10496*a^13*b^4
*c^6 - 14336*a^14*b^2*c^7) + 4*a^11*b*c^8*x)*(-(b^9 - b^4*(-(4*a*c - b^2)^5)^(1/2) + 80*a^4*b*c^4 + 61*a^2*b^5
*c^2 - 120*a^3*b^3*c^3 - a^2*c^2*(-(4*a*c - b^2)^5)^(1/2) - 13*a*b^7*c + 3*a*b^2*c*(-(4*a*c - b^2)^5)^(1/2))/(
512*(a^5*b^8 + 256*a^9*c^4 - 16*a^6*b^6*c + 96*a^7*b^4*c^2 - 256*a^8*b^2*c^3)))^(1/4) + ((-(b^9 - b^4*(-(4*a*c
- b^2)^5)^(1/2) + 80*a^4*b*c^4 + 61*a^2*b^5*c^2 - 120*a^3*b^3*c^3 - a^2*c^2*(-(4*a*c - b^2)^5)^(1/2) - 13*a*b
^7*c + 3*a*b^2*c*(-(4*a*c - b^2)^5)^(1/2))/(512*(a^5*b^8 + 256*a^9*c^4 - 16*a^6*b^6*c + 96*a^7*b^4*c^2 - 256*a
^8*b^2*c^3)))^(3/4)*(4096*a^15*c^8 - x*(-(b^9 - b^4*(-(4*a*c - b^2)^5)^(1/2) + 80*a^4*b*c^4 + 61*a^2*b^5*c^2 -
120*a^3*b^3*c^3 - a^2*c^2*(-(4*a*c - b^2)^5)^(1/2) - 13*a*b^7*c + 3*a*b^2*c*(-(4*a*c - b^2)^5)^(1/2))/(512*(a
^5*b^8 + 256*a^9*c^4 - 16*a^6*b^6*c + 96*a^7*b^4*c^2 - 256*a^8*b^2*c^3)))^(1/4)*(32768*a^16*c^8 + 1024*a^12*b^
8*c^4 - 12288*a^13*b^6*c^5 + 51200*a^14*b^4*c^6 - 81920*a^15*b^2*c^7) + 256*a^11*b^8*c^4 - 2816*a^12*b^6*c^5 +
10496*a^13*b^4*c^6 - 14336*a^14*b^2*c^7) - 4*a^11*b*c^8*x)*(-(b^9 - b^4*(-(4*a*c - b^2)^5)^(1/2) + 80*a^4*b*c
^4 + 61*a^2*b^5*c^2 - 120*a^3*b^3*c^3 - a^2*c^2*(-(4*a*c - b^2)^5)^(1/2) - 13*a*b^7*c + 3*a*b^2*c*(-(4*a*c - b
^2)^5)^(1/2))/(512*(a^5*b^8 + 256*a^9*c^4 - 16*a^6*b^6*c + 96*a^7*b^4*c^2 - 256*a^8*b^2*c^3)))^(1/4)))*(-(b^9
- b^4*(-(4*a*c - b^2)^5)^(1/2) + 80*a^4*b*c^4 + 61*a^2*b^5*c^2 - 120*a^3*b^3*c^3 - a^2*c^2*(-(4*a*c - b^2)^5)^
(1/2) - 13*a*b^7*c + 3*a*b^2*c*(-(4*a*c - b^2)^5)^(1/2))/(512*(a^5*b^8 + 256*a^9*c^4 - 16*a^6*b^6*c + 96*a^7*b
^4*c^2 - 256*a^8*b^2*c^3)))^(1/4)*2i - atan((((-(b^9 + b^4*(-(4*a*c - b^2)^5)^(1/2) + 80*a^4*b*c^4 + 61*a^2*b^
5*c^2 - 120*a^3*b^3*c^3 + a^2*c^2*(-(4*a*c - b^2)^5)^(1/2) - 13*a*b^7*c - 3*a*b^2*c*(-(4*a*c - b^2)^5)^(1/2))/
(512*(a^5*b^8 + 256*a^9*c^4 - 16*a^6*b^6*c + 96*a^7*b^4*c^2 - 256*a^8*b^2*c^3)))^(3/4)*(4096*a^15*c^8 + x*(-(b
^9 + b^4*(-(4*a*c - b^2)^5)^(1/2) + 80*a^4*b*c^4 + 61*a^2*b^5*c^2 - 120*a^3*b^3*c^3 + a^2*c^2*(-(4*a*c - b^2)^
5)^(1/2) - 13*a*b^7*c - 3*a*b^2*c*(-(4*a*c - b^2)^5)^(1/2))/(512*(a^5*b^8 + 256*a^9*c^4 - 16*a^6*b^6*c + 96*a^
7*b^4*c^2 - 256*a^8*b^2*c^3)))^(1/4)*(32768*a^16*c^8 + 1024*a^12*b^8*c^4 - 12288*a^13*b^6*c^5 + 51200*a^14*b^4
*c^6 - 81920*a^15*b^2*c^7) + 256*a^11*b^8*c^4 - 2816*a^12*b^6*c^5 + 10496*a^13*b^4*c^6 - 14336*a^14*b^2*c^7) +
4*a^11*b*c^8*x)*(-(b^9 + b^4*(-(4*a*c - b^2)^5)^(1/2) + 80*a^4*b*c^4 + 61*a^2*b^5*c^2 - 120*a^3*b^3*c^3 + a^2
*c^2*(-(4*a*c - b^2)^5)^(1/2) - 13*a*b^7*c - 3*a*b^2*c*(-(4*a*c - b^2)^5)^(1/2))/(512*(a^5*b^8 + 256*a^9*c^4 -
16*a^6*b^6*c + 96*a^7*b^4*c^2 - 256*a^8*b^2*c^3)))^(1/4)*1i - ((-(b^9 + b^4*(-(4*a*c - b^2)^5)^(1/2) + 80*a^4
*b*c^4 + 61*a^2*b^5*c^2 - 120*a^3*b^3*c^3 + a^2*c^2*(-(4*a*c - b^2)^5)^(1/2) - 13*a*b^7*c - 3*a*b^2*c*(-(4*a*c
- b^2)^5)^(1/2))/(512*(a^5*b^8 + 256*a^9*c^4 - 16*a^6*b^6*c + 96*a^7*b^4*c^2 - 256*a^8*b^2*c^3)))^(3/4)*(4096
*a^15*c^8 - x*(-(b^9 + b^4*(-(4*a*c - b^2)^5)^(1/2) + 80*a^4*b*c^4 + 61*a^2*b^5*c^2 - 120*a^3*b^3*c^3 + a^2*c^
2*(-(4*a*c - b^2)^5)^(1/2) - 13*a*b^7*c - 3*a*b^2*c*(-(4*a*c - b^2)^5)^(1/2))/(512*(a^5*b^8 + 256*a^9*c^4 - 16
*a^6*b^6*c + 96*a^7*b^4*c^2 - 256*a^8*b^2*c^3)))^(1/4)*(32768*a^16*c^8 + 1024*a^12*b^8*c^4 - 12288*a^13*b^6*c^
5 + 51200*a^14*b^4*c^6 - 81920*a^15*b^2*c^7) + 256*a^11*b^8*c^4 - 2816*a^12*b^6*c^5 + 10496*a^13*b^4*c^6 - 143
36*a^14*b^2*c^7) - 4*a^11*b*c^8*x)*(-(b^9 + b^4*(-(4*a*c - b^2)^5)^(1/2) + 80*a^4*b*c^4 + 61*a^2*b^5*c^2 - 120
*a^3*b^3*c^3 + a^2*c^2*(-(4*a*c - b^2)^5)^(1/2) - 13*a*b^7*c - 3*a*b^2*c*(-(4*a*c - b^2)^5)^(1/2))/(512*(a^5*b
^8 + 256*a^9*c^4 - 16*a^6*b^6*c + 96*a^7*b^4*c^2 - 256*a^8*b^2*c^3)))^(1/4)*1i)/(((-(b^9 + b^4*(-(4*a*c - b^2)
^5)^(1/2) + 80*a^4*b*c^4 + 61*a^2*b^5*c^2 - 120*a^3*b^3*c^3 + a^2*c^2*(-(4*a*c - b^2)^5)^(1/2) - 13*a*b^7*c -
3*a*b^2*c*(-(4*a*c - b^2)^5)^(1/2))/(512*(a^5*b^8 + 256*a^9*c^4 - 16*a^6*b^6*c + 96*a^7*b^4*c^2 - 256*a^8*b^2*
c^3)))^(3/4)*(4096*a^15*c^8 + x*(-(b^9 + b^4*(-(4*a*c - b^2)^5)^(1/2) + 80*a^4*b*c^4 + 61*a^2*b^5*c^2 - 120*a^
3*b^3*c^3 + a^2*c^2*(-(4*a*c - b^2)^5)^(1/2) - 13*a*b^7*c - 3*a*b^2*c*(-(4*a*c - b^2)^5)^(1/2))/(512*(a^5*b^8
+ 256*a^9*c^4 - 16*a^6*b^6*c + 96*a^7*b^4*c^2 - 256*a^8*b^2*c^3)))^(1/4)*(32768*a^16*c^8 + 1024*a^12*b^8*c^4 -
12288*a^13*b^6*c^5 + 51200*a^14*b^4*c^6 - 81920*a^15*b^2*c^7) + 256*a^11*b^8*c^4 - 2816*a^12*b^6*c^5 + 10496*
a^13*b^4*c^6 - 14336*a^14*b^2*c^7) + 4*a^11*b*c^8*x)*(-(b^9 + b^4*(-(4*a*c - b^2)^5)^(1/2) + 80*a^4*b*c^4 + 61
*a^2*b^5*c^2 - 120*a^3*b^3*c^3 + a^2*c^2*(-(4*a*c - b^2)^5)^(1/2) - 13*a*b^7*c - 3*a*b^2*c*(-(4*a*c - b^2)^5)^
(1/2))/(512*(a^5*b^8 + 256*a^9*c^4 - 16*a^6*b^6*c + 96*a^7*b^4*c^2 - 256*a^8*b^2*c^3)))^(1/4) + ((-(b^9 + b^4*
(-(4*a*c - b^2)^5)^(1/2) + 80*a^4*b*c^4 + 61*a^2*b^5*c^2 - 120*a^3*b^3*c^3 + a^2*c^2*(-(4*a*c - b^2)^5)^(1/2)
- 13*a*b^7*c - 3*a*b^2*c*(-(4*a*c - b^2)^5)^(1/2))/(512*(a^5*b^8 + 256*a^9*c^4 - 16*a^6*b^6*c + 96*a^7*b^4*c^2
- 256*a^8*b^2*c^3)))^(3/4)*(4096*a^15*c^8 - x*(-(b^9 + b^4*(-(4*a*c - b^2)^5)^(1/2) + 80*a^4*b*c^4 + 61*a^2*b
^5*c^2 - 120*a^3*b^3*c^3 + a^2*c^2*(-(4*a*c - b^2)^5)^(1/2) - 13*a*b^7*c - 3*a*b^2*c*(-(4*a*c - b^2)^5)^(1/2))
/(512*(a^5*b^8 + 256*a^9*c^4 - 16*a^6*b^6*c + 96*a^7*b^4*c^2 - 256*a^8*b^2*c^3)))^(1/4)*(32768*a^16*c^8 + 1024
*a^12*b^8*c^4 - 12288*a^13*b^6*c^5 + 51200*a^14*b^4*c^6 - 81920*a^15*b^2*c^7) + 256*a^11*b^8*c^4 - 2816*a^12*b
^6*c^5 + 10496*a^13*b^4*c^6 - 14336*a^14*b^2*c^7) - 4*a^11*b*c^8*x)*(-(b^9 + b^4*(-(4*a*c - b^2)^5)^(1/2) + 80
*a^4*b*c^4 + 61*a^2*b^5*c^2 - 120*a^3*b^3*c^3 + a^2*c^2*(-(4*a*c - b^2)^5)^(1/2) - 13*a*b^7*c - 3*a*b^2*c*(-(4
*a*c - b^2)^5)^(1/2))/(512*(a^5*b^8 + 256*a^9*c^4 - 16*a^6*b^6*c + 96*a^7*b^4*c^2 - 256*a^8*b^2*c^3)))^(1/4)))
*(-(b^9 + b^4*(-(4*a*c - b^2)^5)^(1/2) + 80*a^4*b*c^4 + 61*a^2*b^5*c^2 - 120*a^3*b^3*c^3 + a^2*c^2*(-(4*a*c -
b^2)^5)^(1/2) - 13*a*b^7*c - 3*a*b^2*c*(-(4*a*c - b^2)^5)^(1/2))/(512*(a^5*b^8 + 256*a^9*c^4 - 16*a^6*b^6*c +
96*a^7*b^4*c^2 - 256*a^8*b^2*c^3)))^(1/4)*2i + 2*atan((((-(b^9 - b^4*(-(4*a*c - b^2)^5)^(1/2) + 80*a^4*b*c^4 +
61*a^2*b^5*c^2 - 120*a^3*b^3*c^3 - a^2*c^2*(-(4*a*c - b^2)^5)^(1/2) - 13*a*b^7*c + 3*a*b^2*c*(-(4*a*c - b^2)^
5)^(1/2))/(512*(a^5*b^8 + 256*a^9*c^4 - 16*a^6*b^6*c + 96*a^7*b^4*c^2 - 256*a^8*b^2*c^3)))^(3/4)*(4096*a^15*c^
8 - x*(-(b^9 - b^4*(-(4*a*c - b^2)^5)^(1/2) + 80*a^4*b*c^4 + 61*a^2*b^5*c^2 - 120*a^3*b^3*c^3 - a^2*c^2*(-(4*a
*c - b^2)^5)^(1/2) - 13*a*b^7*c + 3*a*b^2*c*(-(4*a*c - b^2)^5)^(1/2))/(512*(a^5*b^8 + 256*a^9*c^4 - 16*a^6*b^6
*c + 96*a^7*b^4*c^2 - 256*a^8*b^2*c^3)))^(1/4)*(32768*a^16*c^8 + 1024*a^12*b^8*c^4 - 12288*a^13*b^6*c^5 + 5120
0*a^14*b^4*c^6 - 81920*a^15*b^2*c^7)*1i + 256*a^11*b^8*c^4 - 2816*a^12*b^6*c^5 + 10496*a^13*b^4*c^6 - 14336*a^
14*b^2*c^7)*1i + 4*a^11*b*c^8*x)*(-(b^9 - b^4*(-(4*a*c - b^2)^5)^(1/2) + 80*a^4*b*c^4 + 61*a^2*b^5*c^2 - 120*a
^3*b^3*c^3 - a^2*c^2*(-(4*a*c - b^2)^5)^(1/2) - 13*a*b^7*c + 3*a*b^2*c*(-(4*a*c - b^2)^5)^(1/2))/(512*(a^5*b^8
+ 256*a^9*c^4 - 16*a^6*b^6*c + 96*a^7*b^4*c^2 - 256*a^8*b^2*c^3)))^(1/4) - ((-(b^9 - b^4*(-(4*a*c - b^2)^5)^(
1/2) + 80*a^4*b*c^4 + 61*a^2*b^5*c^2 - 120*a^3*b^3*c^3 - a^2*c^2*(-(4*a*c - b^2)^5)^(1/2) - 13*a*b^7*c + 3*a*b
^2*c*(-(4*a*c - b^2)^5)^(1/2))/(512*(a^5*b^8 + 256*a^9*c^4 - 16*a^6*b^6*c + 96*a^7*b^4*c^2 - 256*a^8*b^2*c^3))
)^(3/4)*(4096*a^15*c^8 + x*(-(b^9 - b^4*(-(4*a*c - b^2)^5)^(1/2) + 80*a^4*b*c^4 + 61*a^2*b^5*c^2 - 120*a^3*b^3
*c^3 - a^2*c^2*(-(4*a*c - b^2)^5)^(1/2) - 13*a*b^7*c + 3*a*b^2*c*(-(4*a*c - b^2)^5)^(1/2))/(512*(a^5*b^8 + 256
*a^9*c^4 - 16*a^6*b^6*c + 96*a^7*b^4*c^2 - 256*a^8*b^2*c^3)))^(1/4)*(32768*a^16*c^8 + 1024*a^12*b^8*c^4 - 1228
8*a^13*b^6*c^5 + 51200*a^14*b^4*c^6 - 81920*a^15*b^2*c^7)*1i + 256*a^11*b^8*c^4 - 2816*a^12*b^6*c^5 + 10496*a^
13*b^4*c^6 - 14336*a^14*b^2*c^7)*1i - 4*a^11*b*c^8*x)*(-(b^9 - b^4*(-(4*a*c - b^2)^5)^(1/2) + 80*a^4*b*c^4 + 6
1*a^2*b^5*c^2 - 120*a^3*b^3*c^3 - a^2*c^2*(-(4*a*c - b^2)^5)^(1/2) - 13*a*b^7*c + 3*a*b^2*c*(-(4*a*c - b^2)^5)
^(1/2))/(512*(a^5*b^8 + 256*a^9*c^4 - 16*a^6*b^6*c + 96*a^7*b^4*c^2 - 256*a^8*b^2*c^3)))^(1/4))/(((-(b^9 - b^4
*(-(4*a*c - b^2)^5)^(1/2) + 80*a^4*b*c^4 + 61*a^2*b^5*c^2 - 120*a^3*b^3*c^3 - a^2*c^2*(-(4*a*c - b^2)^5)^(1/2)
- 13*a*b^7*c + 3*a*b^2*c*(-(4*a*c - b^2)^5)^(1/2))/(512*(a^5*b^8 + 256*a^9*c^4 - 16*a^6*b^6*c + 96*a^7*b^4*c^
2 - 256*a^8*b^2*c^3)))^(3/4)*(4096*a^15*c^8 - x*(-(b^9 - b^4*(-(4*a*c - b^2)^5)^(1/2) + 80*a^4*b*c^4 + 61*a^2*
b^5*c^2 - 120*a^3*b^3*c^3 - a^2*c^2*(-(4*a*c - b^2)^5)^(1/2) - 13*a*b^7*c + 3*a*b^2*c*(-(4*a*c - b^2)^5)^(1/2)
)/(512*(a^5*b^8 + 256*a^9*c^4 - 16*a^6*b^6*c + 96*a^7*b^4*c^2 - 256*a^8*b^2*c^3)))^(1/4)*(32768*a^16*c^8 + 102
4*a^12*b^8*c^4 - 12288*a^13*b^6*c^5 + 51200*a^14*b^4*c^6 - 81920*a^15*b^2*c^7)*1i + 256*a^11*b^8*c^4 - 2816*a^
12*b^6*c^5 + 10496*a^13*b^4*c^6 - 14336*a^14*b^2*c^7)*1i + 4*a^11*b*c^8*x)*(-(b^9 - b^4*(-(4*a*c - b^2)^5)^(1/
2) + 80*a^4*b*c^4 + 61*a^2*b^5*c^2 - 120*a^3*b^3*c^3 - a^2*c^2*(-(4*a*c - b^2)^5)^(1/2) - 13*a*b^7*c + 3*a*b^2
*c*(-(4*a*c - b^2)^5)^(1/2))/(512*(a^5*b^8 + 256*a^9*c^4 - 16*a^6*b^6*c + 96*a^7*b^4*c^2 - 256*a^8*b^2*c^3)))^
(1/4)*1i + ((-(b^9 - b^4*(-(4*a*c - b^2)^5)^(1/2) + 80*a^4*b*c^4 + 61*a^2*b^5*c^2 - 120*a^3*b^3*c^3 - a^2*c^2*
(-(4*a*c - b^2)^5)^(1/2) - 13*a*b^7*c + 3*a*b^2*c*(-(4*a*c - b^2)^5)^(1/2))/(512*(a^5*b^8 + 256*a^9*c^4 - 16*a
^6*b^6*c + 96*a^7*b^4*c^2 - 256*a^8*b^2*c^3)))^(3/4)*(4096*a^15*c^8 + x*(-(b^9 - b^4*(-(4*a*c - b^2)^5)^(1/2)
+ 80*a^4*b*c^4 + 61*a^2*b^5*c^2 - 120*a^3*b^3*c^3 - a^2*c^2*(-(4*a*c - b^2)^5)^(1/2) - 13*a*b^7*c + 3*a*b^2*c*
(-(4*a*c - b^2)^5)^(1/2))/(512*(a^5*b^8 + 256*a^9*c^4 - 16*a^6*b^6*c + 96*a^7*b^4*c^2 - 256*a^8*b^2*c^3)))^(1/
4)*(32768*a^16*c^8 + 1024*a^12*b^8*c^4 - 12288*a^13*b^6*c^5 + 51200*a^14*b^4*c^6 - 81920*a^15*b^2*c^7)*1i + 25
6*a^11*b^8*c^4 - 2816*a^12*b^6*c^5 + 10496*a^13*b^4*c^6 - 14336*a^14*b^2*c^7)*1i - 4*a^11*b*c^8*x)*(-(b^9 - b^
4*(-(4*a*c - b^2)^5)^(1/2) + 80*a^4*b*c^4 + 61*a^2*b^5*c^2 - 120*a^3*b^3*c^3 - a^2*c^2*(-(4*a*c - b^2)^5)^(1/2
) - 13*a*b^7*c + 3*a*b^2*c*(-(4*a*c - b^2)^5)^(1/2))/(512*(a^5*b^8 + 256*a^9*c^4 - 16*a^6*b^6*c + 96*a^7*b^4*c
^2 - 256*a^8*b^2*c^3)))^(1/4)*1i))*(-(b^9 - b^4*(-(4*a*c - b^2)^5)^(1/2) + 80*a^4*b*c^4 + 61*a^2*b^5*c^2 - 120
*a^3*b^3*c^3 - a^2*c^2*(-(4*a*c - b^2)^5)^(1/2) - 13*a*b^7*c + 3*a*b^2*c*(-(4*a*c - b^2)^5)^(1/2))/(512*(a^5*b
^8 + 256*a^9*c^4 - 16*a^6*b^6*c + 96*a^7*b^4*c^2 - 256*a^8*b^2*c^3)))^(1/4) - 1/(a*x)